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Abstract 

Stochastic Loewner Evolution (SLE K ) has been introduced as a description of 
the continuum limit of cluster boundaries in two-dimensional critical systems. 
We show that the problem of ./V radial SLEs in the unit disc is equivalent 
to Dyson's Brownian motion on the boundary of the disc, with parameter 
(3 = 4/k. As a result various equilibrium critical models give realisations 
of circular ensembles with (3 different from the classical values of 1,2 and 
4 which correspond to symmetry classes of random U(iV) matrices. Some 
of the bulk critical exponents are related to the spectrum of the associated 
Calogero-Sutherland hamiltonian. The main result is also checked against 
the predictions of conformal field theory. An erratum to the published 
version has been added in which the above comparison is made 
more explicit. 
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Introduction 



Recently a new method for understanding the scaling limit of conformally invariant two- 
dimensional critical systems has been introduced by Schramm [1] and developed by Lawler, 
Schramm and Werner [2] (LSW). This is known as stochastic Loewner evolution (SLE). It 
relies on the fact that many such systems may be realised geometrically in terms of sets of 
random curves, whose statistics can be described by a stochastic dynamical process. SLE is 
the continuum limit of this process. 

There is in fact a continuous family of SLE processes, labelled by a real parameter k > 0. 
Different values of k are supposed to correspond to different universality classes of critical 
phenomena. For example, for 4 < k < 8 they describe the perimeters of the Fortuin-Kastelyn 
clusters of the Q-state Potts model with 4 > Q > 0, while for 2 < k < 4 they describe the 
graphs of the high-temperature expansion of the O(n) model with — 2 < n < 2 (dual to the 
boundaries of critical Ising spin clusters for n = 1), as well as the external boundaries of 
F-K clusters. This correspondence has so far been proven rigorously only in a few cases [3] . 
However, if the continuum limit of the lattice curves exists and is conformally invariant, it 
must be described by SLE [2]. 

Under these assumptions, LSW [2] have rederived many of the known results for two- 
dimensional critical behaviour which have been found by less rigorous approaches such as 
Coulomb gas methods and conformal field theory, as well as establishing some new ones. 
One aspect of the connection with conformal field theory has recently been pointed out by 
Bauer and Bernard [4] and Friedrich and Werner [5]. 

The particular setting we consider in this note is as follows: consider a critical system 
in a disc of radius R, with a puncture at the origin of radius e, in the limit when R is much 
larger than e and the lattice spacing a. Suppose there are exactly N open curves connecting 
the inner and outer boundaries. In addition, there are no open curves which begin and end 
on the outer boundary. See Fig. 1. For example, these curves could be mutually avoiding 
self-avoiding walks [6], or the external boundaries of percolation clusters (both k — |), or 
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the boundaries of critical Ising spin clusters (k — 3.) (In these last two cases we assume that 
the ensemble is conditioned so as to satisfy the above.) Another example, conjectured to 
correspond to to k — 4, is when the curves are the level lines of a two-dimensional crystalline 
surface at the roughening transition, and there is a screw dislocation of strength N located 
at the origin. As long as k < 4 it is known [7] that, in the continuum limit, these curves are 
simple, that is, they self-intersect with probability zero. For the same reason, the positions 
at which the curves intersect the outer boundary, labelled by complex numbers Re l6 \ are 
well-defined for k < 4. 

Our main result is that, for R ^> e, the joint probability density function (p.d.f.) of these 
points is given by Dyson's circular ensemble [8] 

P eq ({^})oc I] |e^-e^f (1) 

l<j<k<N 

with (3 = A/k. Our argument proceeds by showing that the SLE process appropriate for 
this situation contains the Brownian process invented by Dyson [9], whose equilibrium dis- 
tribution is given by (1), with time being asymptotically proportional to ln(_R/e). 

The distribution (1) is known to describe the statistics of the eigenvalues of random 
unitary matrices in the orthogonal, unitary and symplectic ensembles for (3 = 1,2 and 4 
respectively. Our arguments thus provide simple physical realisations of this ensemble for 
other values of f3 > 1, for example (3 = § (self-avoiding walks) and (3 = | (Ising spin cluster 
boundaries.) 

We also check (1) against the predictions of conformal field theory (CFT). There is a 
subtle factor of \ in the exponent which we elucidate. Dyson's process is known to be related 
by a similarity transformation to the quantum Calogero-Sutherland model [10]. We point 
an interesting connection between the dilatation operator D = L + L of CFT and the 
Calogero-Sutherland hamiltonian. 

It turns out that the eigenvalues of this hamiltonian may, with suitable boundary con- 
ditions, correspond to bulk scaling exponents of the models with k > 4, an example being 
the one-arm exponent computed by LSW [11]. 



Multiple SLEs 



Our arguments are based on an iV-particle generalisation of radial SLE. A single radial 
SLE describes the continuum limit of a curve in the unit disc U := {z : \z\ < 1} which 
begins on the boundary at time t — and ends up at the origin as t — > oo. Let K(t) be 
the hull of the process up to time t (for k < 4 this is just the set of points on the curve.) 
There is a conformal mapping g t : U \ K(t) — > U, such that g t (0) = and g' t (0) > 0. LSW 
argue that gt(z) may be chosen so as to satisfy an evolution equation 

g t (z) + e^ B ^ 

g t {z) = -g t (z) gt(z) _ ^ KB{t) (2) 

where B(t) is a standard one-dimensional Brownian process with E[£?(t) 2 ] = t. Note that 
time has been reparametrised so that <^(0) = e*. (2) is the standard form of radial SLE, 
which maps the trace of the SLE into the point e l V*-B{t) on ^ e boundary, but for our purposes 
it is more convenient to consider g t (z;9) = g t (z)e^ e ~^ B ^\ which maps the trace into e 10 , 
and satisfies 1 

dg t (z; 6) = -g t {z- 6) dt - ig t (z; 9)^dB{t) (3) 

Now consider N SLE's which start from distinct points {e l6lj } on the boundary, with 
1 < j < N. Let Kj{t) be the hull of the jth SLE. For k < 4 these are segments of non- 
intersecting simple curves. Let G[ N \z) be a function which conformally maps U\uf =1 Kj(t) 
onto U, with ^(0) = and G^^O) > 0. Then we shall argue that G { t N \z) may be 



chosen to satisfy 



G\ N) = -G[ N) £ G * +e (4) 

^ _ e iOj(t) 



where 



1 In the Ito convention there is an additional term —^gdt on the right hand side. This disappears 
again after making the global rotation leading from (8) to (4). 
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dBj(t) = cot((^(t) - 9 k {t))/2) dt + y/KdBj(t) (5) 

and Bj(t) are N independent Brownian motions, starting at the origin. 

To see this, consider the infinitesimal transformation o (G^)^ 1 and note that this 
may be obtained by allowing each SLE to evolve independently according to (3) over a time 
dt: 

G[ N + i o (G?V = g£\9 N (t)) o g^iO^it)) o . . . o g^M)) (6) 

During the evolution of the jth SLE, G[ N ^ evolves according to (3), with 6 = 9j(t) and 
B(t) = Bj(t), but the other 9 k (t) with k ^ j also evolve according to 

e iO k (t) i e i0j(t) 

d6k ^ = 1 e ie k (t) _ e ie 3 (t) dt - y/KdBj(t) = cot((6 k (t) - 9 3 {t))/2)dt - yfcdBj(t) (7) 
Thus, after evolving every SLE with j — 1, . . . , N, we have 

N r {N) i6 u) N 



dC'r = -GT } Y dt - iGT^YdBAt) (8) 



where 



d0,-(t) = ^cot((^(t) - 9 k (t))/2)dt - ^Y. dMt) (9) 

It is now simpler to rotate the whole disc through an angle v^SjLi dBj(t), after which we 
obtain (4) and (5) as claimed. 

Eq. (5) is the Dyson process [9]. It may be written as 

dV ,- 

d6j = -—dt + y/KdBj(t) (10) 

where V = — 2X) J - <fc ln | sin((0j — 9 k )/2)\. At late times, the distribution of the {9j(t)} tends 
towards an equilibrium p.d.f. at temperature n/2: 

P cq (9 u ...,6 N )cx e- 2V ' K = I sin((^ - 9 k )/2)\ A ' K oc J] \e»* - e^f (11) 

j<k j<k 

This is Dyson's circular ensemble [8], with 
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/3 = 4/« . (12) 

So far, the starting points {#,} have been taken to be fixed. Now consider an ensemble 
of these, generated by the configurations of some bulk critical ensemble in the interior of 
U (conditioned if necessary on the existence of exactly N curves.) We argue that if this 
corresponds to a conformally invariant bulk system, the p.d.f. of the {9j} must be given by 
P cq . Let K({6j},t) be the union of the N hulls up to time t, given that they start at the 
points {e j }. The expectation value of any observable may be taken by first conditioning on 
the subset K({9j},t'), with t' < t. By conformal invariance, the distribution of K({9j},t) \ 
K({8j}, t') is the same as that of its image under G^\ namely K({8j(t')}, t — t'). Averaging 
over K({9j},t') is equivalent to averaging over the {Bj} up to time t' . Taking t — > oo, we 
conclude that 



P{{9 3 }) = E {Bj{t „y, //mt , ]} [P({0,-(f )})] (13) 



that is P({6j}) is stationary under the process (5), and must therefore be equal to P cq . 

This is strictly valid only when the SLEs are allowed to reach the origin. To discuss the 
case when they reach only the circle \z\ — e, it is helpful to map conformally the annulus 
to a cylinder of length I = ln(i?/e). The points e j6lj are now arrayed around one end of the 
cylinder. Since G[ N) = NG { t N \l + 0(Gf°)) as G[ N) -> 0, we see that as long as I > 1 the 
effect of the evolution is to reduce the length of the cylinder at a rate t = —N. Meanwhile the 
points {9j} are moving according to (5). The approach of their distribution to equilibrium 
is expected to be exponentially fast with a rate constant 0(1). Thus, as long as £ 3> N, 
we may apply the same argument as above, and deduce that the distribution of the {9j} is 
given by P cq , with corrections suppressed by powers of e/R. The same should apply on a 
lattice, as long as the spacing a < e <^ R. 
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Comparison with conformal field theory 

The crucial assumption of conformal invariance made in deriving the above result would 
appear to be stronger than the analogous statement for AT = 1. In particular, it is not clear 
from this point of view why invariance under the uniformising transformation G\ N \ which 
assumes that the curves grow at the same rate, is to be chosen among other possibilities, 
although it appears to be the most natural one. For this reason we have checked our main 
result (1) using methods of CFT. In this language we expect the joint p.d.f. to be given by 
the correlation function in the 0{n) conformal field theory 

($ a ,.. OjV (O)0 ai (it:e^) . . .<j> aN (Re ieN )) (14) 

where 4> a (Re ld ) is a boundary 1-leg operator carrying 0(n) index a, and $ is a bulk iV-leg 
operator. By choosing the aj to be all different, we ensure that the curves all reach the 
origin without annihilating. In [12] it was conjectured that the operators (j) a correspond 
to Virasoro representations labelled by (1,2) in the Kac classification. These have a null 
state at level 2, and therefore their correlators satisfy second-order linear partial differential 
equations with respect to each of the 0j (the BPZ equations [13]). The general solution for 
such an iV + 1-point correlator is not however known. Instead, we may take the form (1) as 
an ansatz, and check whether it satisfies these equations. Even this is somewhat tedious, and 
we have carried it through only for iV = 2. Alternatively, one may check whether (1) satisfies 
the fusion rules which follow from the BPZ equations. These determine the behaviour of 
the correlator (14) in the limits when (say) p of the 9j approach each other. Suppose, for 
example, that \9j — 9 k \ = 0(5) for 1 < j < p and 1 < k < p, with 2 < p < N. In the limit 
5 — > we may use the operator product expansion (OPE) 

n0 aj (^)oc^-^0 ai ... ap (^) (is) 

where <t> ai ...a p is the boundary p-leg operator, and x p is its scaling dimension. Given that 
the 1-leg operator corresponds to (1,2), the fusion rules determine the allowed values of 
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Xp which may occur on the right hand side of (15): they are the scaling dimensions /ii,n+i 
of the + 1) operators, with < n < p and p — n even. Duplantier and Saleur [14] 
argued that the p-leg operator must in fact correspond to n — p. Using the Kac formula 
hi, p +i = p{2p + 4 — k)/2k then gives the exponent in (15) to be simply p{p — 1)/k. 

On the other hand, we may simply take the appropriate limit in the ansatz (1) to find 
a dependence on 5 of the form ni<j<fc<p<^ — 5 p ( p ~ 1 ) /3 / 2 . Comparing these two expressions 
apparently gives (5 — 2/k, not A/k as found above. One may take further short-distance 
limits within (1): the results all consistent with Duplantier and Saleur [14] but only if 
(3 = 2/k 2 

The resolution of this paradox is as follows: the correlator (1) may be written in operator 
language as 

(®N\e- w \{9 3 }) (16) 

where t = ln(i?/e), D = L + L is the generator of scale transformations, and \{0j}) is a 
boundary state. The usual formalism of CFT assumes that D is self-adjoint: this originates 
in the invariance of the bulk theory under inversions z — > —1/z. Equivalently, on the 
cylinder, e~ aD is the continuum limit of the transfer matrix, which, for many lattice models, 
may also be chosen to be self-adjoint. But for the loop representation of the 0(n) model 
this is not the case: the only practicable transfer matrix which has been employed [15] acts, 
at a given time t, on the space spanned by a basis defined by the positions of the points at 
which the loops intersect the chosen time slice, together with their relative connections in 
the 'past', but not the 'future'. This asymmetry leads to a transfer matrix T which is not 



2 We remark that the reason that the iV + 1-point correlator (14) has such a simple form is that 
there is only ever a single term on the right-hand side of OPEs like (15). This is because all the 
0(n) indices on the left-hand side are different, and therefore the operators on the right hand side 
transform according to the totally symmetric representation of 0(n). The other possible terms, 
with n < p, correspond to other representations, and could only arise if some of the cij were equal. 
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self-adjoint. In our case, the points at which the loops intersect a given time slice can either 
be connected back to the points e ldj at t = 0, or to each other via loops which close in the 
past at times t > (see Fig. 2). Let IT be the projection operator in the full space in which 
T acts, onto the subspace spanned by the boundary states | {#;}). That is, II traces over 
the positions of the points not connected to t — 0. Then T = HT acts wholly within this 
subspace. Note that JUjd9j({9j}\ is a left eigenstate of T with unit eigenvalue. 
As usual, we may also think of T as acting on L 2 functions of the {9j} through 

(Tfm}) = I udom}\ f m})m}) ( i7 ) 
j j 

so that T^l = 1. Moreover, as ki(R/e) — > oo, the joint p.d.f. of the {9j} is given by the 
right eigenfunction of T with the largest eigenvalue. 

On the other hand, the Langevin equation (5) yields the Fokker-Planck equation P = CP 
for the evolution of P({9j},t) where 

N d dV k d 2 
C= ^ 1 W j W J + 2d9 2 (18) 

where CP eq = and C^l = 0. Now the conformal mapping G[ N ^ oc e Nt acts as a scale 
transformation near the origin, while its action on the unit disc is described by C We 
therefore conjecture that this is nothing but the continuous version of T, more precisely, 
ITT* ~ e atc/N ^ ac tj n g on foe subspace. 

Now for any dynamics which satisfies detailed balance C is related to a self-adjoint 
operator H by a similarity transformation: 

H = -P-V*CPW (19) 

The ground state eigenfunction of H is then P^ 2 . This square root is the origin of the 
discrepancy between the CFT result and (1) with f3 — 4//c: the self-adjoint operator e~ Ht is 
proportional to the scale transformation operator e~ Dt where D = L + L of CFT, acting 
on the subspace spanned by the boundary states, that is e~ m l N oc He -D 'n''\ But it does 
not give the continuum limit of the correct transfer matrix. By assuming that this limit was 
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self-adjoint, which is assumed in the standard formulation of CFT, we found, erroneously, 
the square root of the correct result (1) with (3 = 4/ ' k. 

In the case of Dyson's brownian motion the hamiltonian H is that of the quantum 
Calogero-Sutherland model [10] 

d 2 | 2-« v 1 N(N-l) 

2 ^ 89] + 2k f£ <** 2 (0j ~ 0k)/2 2k { U) 

It is interesting to note that the adjoint operator C) = e 2V ^ K Ce~ 2V ^ K has the form 

£t = f E^ + EE ™m - e k )/2)± (2i) 

3 3 j k^j 3 

C) is the generator for a typical first-passage problem. For example, in the case N = 2, 
the probability h(9i, 9 2 ; t) that the two particles have not met up to time t, given that they 
started from (9i,9 2 ) satisfies 8 t h = C)h. In fact, with 9 = Q\ — 9 2 , and rescaling 2t — > t, 
this is just the equation derived by LSW [11], whose lowest non-trivial eigenvalue gives the 
one-arm exponent, related to the fractal dimension of F-K clusters for k = 6. Eigenfunctions 
of C) behave near 9 = as 9 a where a = or 1 — 4/ k. When k < 4 the appropriate solution 
corresponds to a — 0, or h — 1, consistent with the result that SLE is a simple curve, 
but when k > 4 the solution is non-trivial. LSW [11] argue that the appropriate boundary 
condition at 9 = 2-7T for the one-arm problem is dh/89 = 0, and that the solution is then 
( sin(^/4)) 1 ~ 4 ^ K e~ xt with A = (k 2 — 16)/32k. This is of course also an eigenvalue of the 
Calogero-Sutherland hamiltonian (20), and it raises the question as to whether other bulk 
scaling dimensions are given by eigenvalues of Calogero-Sutherland systems with suitable 
boundary conditions. 
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Erratum to the published version. 

In the published paper [16] (as its appears above), it was noted that the only form of 
the result for the joint p.d.f. of the boundary points that is consistent with conformal field 
theory (CFT) is Y\j < k(e t9j — e dk ) 2 / K , which is different from the equilibrium distribution of 
the corresponding Dyson process, where the exponent is 4/k. A possible explanation of this 
discrepancy was given. 

Recently [17], however, we have performed an ab initio CFT calculation. This confirms 
the exponent 2/k, but it also shows the correct source of the discrepancy lies in the assump- 
tion above Eq. (13) that the measure on the curves is conformally invariant. This is too 
strong - if instead we allow for it to be invariant up to a conformal factor n? \g't( etej )\ h2,1 
(where /i2,i = (6 — k)/2k), the results of the two computations agree. 

The corresponding Calogero-Sutherland model then turns out out to have (5 — 8/k. The 
reader is referred to [17] for details. 
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FIGURES 




FIG. 1. The geometrical set-up. N (here = 3) open curves connect the boundaries at r = e 
and R of the annulus, intersecting the outer boundary at points {Re 1 ® 3 }. No other open curves 
are allowed to intersect r = R, but they may intersect r = e, as well as there being any number of 
closed loops (except when they carry zero weight, as for n = 0.) 




^ t 

FIG. 2. The transfer matrix of Ref. 15 on the cylinder. It keeps track of the positions of points 

where curves intersect the given time-slice, as well as their connectivity in the past (but not the 
future.) We distinguish between those points connected to the boundary at t = (r = R), and 
those which close for times t > 0. 
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